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1 Introduction 

Trace formulae and the non-linear supersymmetric cx-model are basic analytical 
tools used successfully in the fields of quantum chaos and disordered systems. Both 
are designed to treat systems with a small number of degrees of freedom. Hence they 
are limited in their possibility of analyzing many-body systems where interparticle 
interactions play an important role, and the number of degrees of freedom is large. 

On the other hand, many experimental studies of quantum chaos use systems 
which consist in a large number of interacting particles, for example quantum dots or 
disordered metallic particles. Having an elaborate single-particle description of these 
systems, it is of prominent importance to understand the role of interactions, the range 
of applicability of a single-particle picture, and the interplay between chaos and inter- 
particle interactions. 

In this respect, an important observation is that strong chaotic dynamics, on the 
level of non-interacting single-particle description, provides us with the possibility of 
analyzing interacting many-body systems by a systematic perturbative approach. The 
small parameter of this perturbation theory is 1/g, where g = tn/t c is the dimensionless 
conductance, i.e. the ratio of the Heisenberg time, t# (the inverse mean level spacing), 
to the classical relaxation time, t c . 

The general form of the interaction Hamiltonian in which particles interact via a 
two-body potential U(r, r') is 



where c ia and q ct are the creation and annihilation operators for a particle in state ipi 
and spin a, while 



. _ 1 
int o 




ijkl era' 




I 



drdr'U{r, r') V* (r) ip* (r')^ (r) ipk (r') 



are the matrix elements of the interaction potential. These matrix elements can be 
divided into two groups according to their typical magnitude. One contains diagonal 
matrix elements, namely those Uijki in which two pairs of indices are identical. All the 
other matrix elements, which we call off-diagonal, are included in the second group. In 
appendix A it is shown that the typical magnitude of off-diagonal matrix elements is as 
small as d/g where d is the single-particle mean level spacing and g is the dimensionless 
conductance. The same smallness restricts also the fluctuations in the diagonal matrix 
elements. Therefore, the interaction Hamiltonian of electrons in a quantum dot takes 
the form 
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The first term of this formula is the orthodox model 1} 2 representing the charging energy 
of the dot: C is the total capacitance of the dot, iV = c la c ia- is the total number of 
electrons in the dot, and N is a continuous parameter controlled by the gate voltage. 
The second term, when A > 0, represents an attractive interaction which drives the 
grain into a superconducting state at sufficiently low temperatures and weak magnetic 
field. The last term represents the electron-electron interaction in the spin channel, 
and the coefficient a is of order the single-particle mean level spacing, d. 

Strong chaotic dynamics of the noninteracting particles implies that 
g ^> 1. Consequently, all off-diagonal matrix elements Uiju are proportional to 1/g, 
and the mean field approximation (|2p for the interacting Hamiltonian (|I]) is justified. 
Indeed, many phenomena of normal and superconducting metallic grains are described 
by the approximation (0). The most prominent one is the Coulomb blockade, which 
is essentially the quantization of the number of electrons in the grain away from the 
charge degeneracy point. Because of this quantization, the zero-bias conductance of the 
system vanishes, while the current / as the function of the source-drain voltage V shows 
a threshold behavior. The fine structure of the current-voltage curve is associated with 
the single-electron levels of the system 3 ' 4 . 

Nevertheless, there are interesting phenomena emerging from the fluctuations of 
the interaction matrix elements, i.e. with the 0(d/g) corrections to (0). In this review 
we analyze two experiments of Ralph, Black and Tinkham 5 ' 6 and show that fluctuations 
of the interaction energy, although small as d/g, clearly manifest themselves in the 
differential conductance spectra of ultrasmall metallic grains. The small effects of 
fluctuations in the charging energy are especially pronounced due to the fact that the 
system is driven out of equilibrium, and is able to explore several high excited states 
at relatively low source-drain voltage. 

The experimental system consists of a single aluminum particle connected to ex- 
ternal leads via high resistance (1-5 MO) tunnel junctions formed by oxidizing the 
surface of the particle. The device (see illustration in Fig. 1) is fabricated using electron 
beam lithography and reactive ion etching to form a bowl-shaped hole in an insulating 
S13N4 membrane. The opening at the lower edge of the membrane is 3-10 nm in diam- 
eter. Al is evaporated onto the bowl-shaped side of the membrane, and subsequently 
the Al surface is oxidized. The oxide layer forms a tunnel barrier in the vicinity of 
the small hole in the membrane. The membrane is then flipped up side down and a 
small amount of Al is deposited. Because of surface tension, the Al forms a layer of 
electrically isolated particles, a few nanometers in size. Following a second oxidation, 
a thick layer of Al is deposited on top of the particles.^ In approximately 25% of the 

* In other configurations of these devices a gate electrode of ring shape is deposited after flipping the 
membrane, and the same procedure follows the oxidation of the gate. 
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Figure 1. A schematic illustration of the device used in the experiments of Ralph, Black and 
Tinkham for measuring the differential conductance spectra of ultrasmall aluminum grains. 

devices one Al particle covers the hole in the nitride membrane, so that the electrons 
passing between the leads tunnel trough the metal particle. 

The capacitances and the resistances of the tunnel junctions are estimated by 
fitting the large scale I — V curves, eV ~ e 2 /C, to the Coulomb blockade staircase 
pattern. From the capacitances one can determine the area of the tunnel junctions and 
the volume of the grain which is used in turn to estimate the single-particle mean level 
spacing 5 . 

In this review we focus our attention on scales of the I — V curves which are much 
smaller than that of the Coulomb blockade, scales over which the single-particle mean 
level spacing, d, and the fluctuation in the charging energy, d/g, are resolved. Fig. 2 
displays the differential conductance, dl/dV, of two different normal metallic particles 
(of sizes roughly 2.5 and 4.5 nm) as a function of the source-drain bias energy eV. The 
spectra display three clear features: 

1. The low resonances of the differential conductance are grouped in clusters. The 
distance between nearby clusters is of order the mean level spacing d of the 
noninteracting electrons in the dot. 

2. The first cluster contains only a single resonance. 

3. Higher clusters consist of several resonances spaced much more closely than d. 

In section 2 it will be shown that these features are manifestations of the interplay 
between electron-electron interactions and nonequilibrium effects 7 . Each cluster of res- 
onances is identified with one excited single-electron state, and each resonance in turn is 
associated with a different occupancy configuration of the metal particle's other single- 
electron states. The appearance of multiple resonances reflects the strongly nonequi- 
librium state of the particle. 

In another experiment 6 , Ralph, Black and Tinkham measured the tunneling res- 
onance spectra of ultrasmall superconducting grains. The number of electrons in the 
system was controlled by a gate voltage. The results of this experiment, depicted in 
Fig. 3, show that: 

1. For the ground state of the grain with an even number of electrons, the first peak 
of the differential conductance is merely shifted by the gate voltage V g . The shape 
of this peak does not change over a large interval of V g . Contrarily, if the grain 
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Figure 2. The low temperature (30 mK) differential conductance dl/dV versus bias energy of 
ultrasmall Al particles with volumes w 40 nm 3 (upper panel) w 100 nm 3 (lower panel). The first 
resonance is isolated while subsequent resonances are clustered in groups. The distance between 
nearby groups of resonances is approximately the single-particle mean level spacing d. (From 
Rcf. [5]). 
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Figure 3. The tunneling resonances of superconducting grains in the odd (upper scans) and the 
even (lower scans) charging states. Different scans correspond to different value of the gate voltage, 
and are artificially shifted in energy to align peaks due to the same eigenstate. In contrast with 
Fig. 2 the first resonance, in the odd charging state, develops a substructure when shifted by the gate 
voltage. (From Ref. [6]). 
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contains an odd number of electrons, the height of the first peak rapidly reduces 
with a change of the gate voltage, and a structure of subresonances develops on 
the low-voltage shoulder of this peak. 

2. The characteristic energy scale between subresonances of the first peak is of the 
order of the mean level spacing d. 

These observations contrast the results for the normal case in which the first peak 
did not split. Nevertheless, it was suggested in Ref. [6] that the substructure of the 
first peak is still associated with nonequilibrium steady states of the grain. In section 
3, the origin of these nonequilibrium states and the mechanism which generates them 
will be clarified 8 . 

The explanations for both experiments discussed here rely on the assumption that 
the systems are stimulated into steady states which are far from equilibrium, namely 
that relaxation processes are too slow to maintain the system in equilibrium. In section 
4 we summarize the various relaxation processes and estimate the inelastic time, Tj n , 
for electrons in the dot. The results will be summarized in section 5. 

2 Normal grains 

Our model for the experimental system is given by the Hamiltonian: H = Ho + 
Ht + Here Ho describes the noninteracting electrons in the left (L) and right (R) 

leads and in the metallic grainf], 

H o= Y.^ q di q d aq + Y.^ c W ( 3 ) 

a=L,R q I 

Tunneling across the barriers is described by 

H t= E E^ a) <Q+ H.c, (4) 

a=L,R q,l 

where T^f* are the tunneling matrix elements. Interaction effects given by (HD are 
taken into account only for the electrons in the grain, but including screening by image 
charges in the leads. 

For the ultrasmall aluminum grains considered here, one can neglect superconduct- 
ing pairing since the single-particle mean level spacing, « 1 meV, is larger than the BCS 
superconducting gap which is 0.18 meV 9 . Under this condition, the interaction term of 
the electrons is generally approximated by the orthodox model, H- in ^ (eJ2i C ] c i) 2 /C, 
where C is the effective capacitance of the grain 1 . Within this approximation the 
charging energy depends only on the total number of electrons in the dot, but not 
on their particular occupancy configuration. The orthodox model is able to account 
for the Coulomb blockade 1 , and the Coulomb staircase behavior of the current as the 
number of extra tunneling electrons in the dot increases. It can also be generalized to 
describe features on the scale of the single-particle level spacing 10 . However, the or- 
thodox model cannot account for the clusters of resonances in Fig. 2, since these result 
from fluctuations, SU, in the interaction energy between pairs of electrons. 

We focus our attention on the (experimental) voltage regime where there is no 
more than one extra tunneling electron in the dot. At small voltage bias, V, within the 
Coulomb-blockade regime [Fig. 4(a)], current does not flow through the system. Current 

^Unless explicitly written, from now on single-particle and spin states will be denoted by a single 
subscript. 
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Figure 4. An illustration of transport through the metal particle at various values of the 
source-drain voltage V. Filled single-particle levels are indicated by full circles and empty ones by 
open circles. U is the charging energy, and d is the single-particle mean level spacing, (a) The system 
at small voltage bias within the Coulomb blockade regime; (b) V corresponding to the first resonance 
in Fig. 2. The thin dashed lines indicate the energy of a level after an electron has tunneled into the 
dot; (c) V near the first cluster of resonances in Fig. 2. The splitting within the first cluster originates 
from the sensitivity of level i + 1 to the different possible occupation configurations as shown. 



first starts to flow when one state % inside the grain becomes available for tunneling 
through the left barrier, say, as illustrated in Fig. 4(b). As the system becomes charged 
with an additional electron, the potential energy of the other electrons in the dot 
increases by U ~ e 2 /C, and some of the lower energy occupied electronic states are 
raised above the right lead chemical potential [in Fig. 4(b) these "ghost" states are 
shown as dashed lines]. Electrons can tunnel out from these states into the right lead 
leaving the particle in an excited state. There is, however, only one configuration of the 
electrons which allows an electron to tunnel into level % from the left lead, namely all 
lower energy levels occupied. This implies that only a single resonance peak appears in 
the differential conductance at the onset of the current flow through the system (broken 
spin degeneracy would cause splitting of this peak). 

The situation changes when V increases such that electrons can tunnel from the 
left lead into the next higher available state % + 1, as shown in Fig.4 (c). In this case, 
there are several possible occupancy configurations, on which the exact energy of level 
i + 1 depends. The several possible energies of level i+1 lead to a cluster of resonances in 
the differential conductance of the grain. The scenario described above holds provided 
that inelastic processes are too slow to maintain equilibrium in the particle. 

To explicitly demonstrate the splitting of resonances induced by fixed fluctuations 
in the interaction energy 5U, model detailed-balance equations 10 were solved numeri- 
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Figure 5. Model differential conductance obtained from noncquilibrium detailed-balance equations: 
solid line - in the absence of inelastic processes, I/t^ = 0; dashed line - with inelastic relaxation 
rate larger than the tunneling rate, 1/Tj n = 5/r t , m . 



cally and the corresponding differential conductance plotted in Fig. 5 by the solid line. 
The model system consists of 7 equally spaced levels, occupied alternately by 4 or 5 
electrons, in a current-carrying steady state. For simplicity, the tunneling rate into 
each level, l/r tun (Tl^(€i) in the notation of Ref. [10]), is chosen to be uniform, and 
the voltage is applied by increasing the left chemical potential. The temperature is 
1% of the mean level spacing d, and the variance of the fluctuations 5U in the inter- 
action energy is d/5. In the absence of fluctuations (SU = 0), dl/dV consists of single 
resonances spaced by d. 

To estimate the fluctuations in the interaction energy consider the Hartree term 
of the interaction energy, Uh- We wish to calculate the interaction energy difference 
associated with different occupation configurations of low energy states. Suppose that, 
as illustrated in Fig. 4(c), these differ by a single occupation number, namely, in one 
configuration the state j is empty and f is full while in the other f is empty and j is 
full. Then 



SUh= /dridr2|Vi(rx)rt7(ri,r a ) |^"(r 2 )| 2 - |^(r 2 



where the index i labels an electron state other than j or j', f/(r 1; r 2 ) is the interaction 
potential. Since wave functions of chaotic systems associated with different energies are 
statistically independent, (5Uh) = where (• • •) denotes ensemble or energy averaging. 
We are therefore interested in fluctuations of SUh which emerge from the non-uniform 
probability distributions of the single-particle eigenstates in real space. The calculation 
of (SUfj) is similar to that presented in appendix A. The result for diffusive systems is 



'9, 



where c = \/2aY, n l n l ~ 4 / 7r is a constant of order unity, and a equals two for system 

with time reversal symmetry and unity for systems without time reversal symme- 
tryQ. The above estimate for the fluctuations in the charging energy also applies for 



■'■This estimation does not take into account a change in the potential due to the insertion of an 
additional electron. It was argued that the latter effect may lead to an even stronger effect, i.e. 
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general chaotic systems, with g = h'ji/d where 71 is the first non- vanishing Perron- 
Frobenius eigenvalue 12 , see appendix A. 

The increase of the fluctuations in the interaction energy as g decreases is related 
to the fact that g is a measure for the uniformity of the single-particle wave functions. 
The bigger g the more uniform are the wave functions and the less are the fluctuations 
in the interaction energy. Experimentally we find g ~ 5. Unfortunately, an analytical 
estimate of g requires precise knowledge of the shape and disorder of the particle which 
we lack. A naive estimate of g in ballistic systems is h/rd, where r is the time for an 
electron at the Fermi energy to cross the system. The metallic grains of the experiment, 
however, have a roughly pancake shape. Assuming diffusive dynamics one can show 
that grd/h oc (z/r) 2 where z is the pancake thickness and r is its radius, g is therefore 
much smaller than h/rd. 

When M available states below the highest accessible energy level (including spin), 
are occupied by M' < M electrons, there are (mO different occupancy configurations. 
The typical width of a cluster of resonances in this case is W 1 l 2 cdjg where W = 
min(M — M', M'). The width of a cluster of resonances therefore increases with the 
source-drain voltage. The distance between nearby peaks of the cluster, on the other 
hand, decreases as W 1 / 2 /(m')- This behavior can be seen in Fig. 5. 

3 Superconducting grains 

As illustrated by Fig. 4, the splitting of tunneling resonance peaks in normal 
metallic grains comes from the possibility of forming different occupation configurations 
of single-particle states at sufficiently high source-drain voltage. These configurations 
are reached by resonant tunneling provided relaxation processes are sufficiently slow. 
This picture also explains the observation that the first cluster in Fig. 2 contains only 
a single peak. However, the data of Fig. 3 shows that the first peak in the tunneling 
resonance spectra also splits into several subresonances (see illustration in Fig. 6). This 
behavior appears when the superconducting grain contains an odd number of electrons 
and the gate voltage is such that the dot is far from the charge degeneracy point. 

In this section we show that the development of a substructure in the first peak of 
the tunneling resonance spectra is also associated with the generation of nonequilibrium 
steady state. However in contrast with the resonant tunneling mechanism used in the 
previous section for the high resonance peaks, here the nonequilibrium steady state is 
reached by inelastic cotunneling processes. 

The principal difference between odd and even grains is that all excitations of the 
latter are of energy larger than the superconducting gap 2A. Therefore, a source-drain 
voltage in the range V < A/e can not induce nonequilibrium states. Odd grains, on 
the other hand, contain one unpaired electron, which may be shifted to various single- 
electron levels with characteristic energy scale smaller than the mean level spacing d. 
For this reason even a small source-drain voltage d < eV < A is sufficient to excite 
the grain. The mechanism of excitation is inelastic cotunneling 13 . Tunneling into the 
excited grain requires less energetic electrons, and lead in turn to the substructure on 
the low- voltage shoulder of the of the first resonance, see Figs. 3 and 6. A closely related 
problem was considered by Averin and Nazarov 14 , however, their theory assumed that 
relaxation processes prevent the formation of nonequilibrium states. As will be argued 
in the next section, relaxation processes in ultrasmall metallic grains are very slow, and 
therefore will be neglected in our theory. 

To describe the effect quantitatively, we construct the master equations governing 
the time evolution of probabilities of different electronic configurations of supercon- 
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Figure 6. A schematic illustration of the differential conductance of an "odd" superconducting 
grain as function of the source-drain voltage V, at various gate voltages V g . Higher resonances are 
separated by the superconducting gap from the first one, and subresonances are developed as the 
first resonance is shifted by the gate voltage. 

ducting grains allowing for second order cotunneling processes. The solution of these 
equations for two limiting cases (one in which two levels participate in the transport, 
and the other when a large number of levels contribute) explains the substructure of 
the first peak of the differential conductance illustrated in Fig. 6. 

As in the previous section the model Hamiltonian consists of three terms: H = 
H + Ht + -f^int- -^Oj given by fl3|), describes the noninteracting electrons in the leads 
and in the dot; Ht, given by (f|), is the tunneling Hamiltonian, and the interaction 
Hamiltonian will be approximated by 



where N = c }a c jcr is the number of electrons in the dot, and N Q is a continuous 
parameter controlled by the gate voltage. Nq determines the finite charging energy 
required to insert, U+, or to remove, LL, one electron, 



Consider the experimentally relevant case, e 2 /C ^> A, so that the grain has well 
defined number of electrons. If this number is even iV = 2m, the ground state energy 
(we will omit charging part of the energy and restore it later), E 2 m = Fi m + 2/zm, 
can be calculated in the mean field approximation 15, 9 by minimizing thermodynamic 
potential F 2m = J2k(£k — £fc) + A 2 /A where = (£j + A 2 ) 1 / 2 , with respect to A, and 
by fixing the chemical potential according to the number of electrons in the grain. All 
the excited states of even dots are separated from the ground state by a large energy, 
2 A. Considering now the energy spectrum of an odd grain, iV = 2m — 1, we notice 
that the second term in Eq. (^j) operates only within spin singlet states. Therefore, to 
calculate the low-lying excited states in this case, we fill the single-electron state j with 
one electron, and then find the ground state of the remaining 2m electrons with state j 
excluded from the Hilbert space. In the mean field approximation it corresponds to the 
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minimization of the thermodynamic potential i^m-i = Hk^ji^k — £fc) + A 2 / A + £j- The 
excited states with energies smaller than A are characterized by a single index, j and 
will be denoted by E^-i- in what follows we will need the energy cost of introducing 
an additional electron into the odd state: U + +Ej, where Ej = -Esm - -E^m-i- 111 appendix 
B it is shown that in the limit A>d the result is: 

3d ^jd 



^ = ^m-- + f K -^e 3 +^. (7) 

We turn now to the kinetics of a superconducting grain. Consider the regime where 

2 2 

U + = U < A, C/_ ~ ^ ^> [/, and ^ 3> A ^> d. We also assume the conductance of the 
tunnel barriers to be much smaller than e 2 /h, and that the source-drain voltage is small 
eV < A. The simplicity brought to the problem in this regime of parameters stems 
from the fact that there is only one available state with an even number of electrons 
(because U- 3> U+ one can only add an electron to grain but not subtract one), and 
whenever the grain contains an even number of electrons it is in its ground state. This 
imply that even grains cannot be driven out of equilibrium state, while for odd grains 
tunneling (and cotunneling) takes place via unique state. 

Henceforth, we concentrate on grains with an odd ground state. Let us denote by 
P e the probability of finding the grain with an even number of electrons, and by Pj the 
probability to find the grain in the odd state j. Since these states are spin degenerate in 
the absence of magnetic field, Pj will denote the sum Pj^ +Pj ±. The master equations 
for the probabilities P e , Pj have the form 

= E [^iePj ~ ^loPe] , (8) 

dP 

_i=2Vrr- p-t -p-i+2r (i) p-r (i) p 

where T^X e an d ^eXo are ^he transitions rates from the odd j-th state to the even and 
from the even to odd respectively, while is the rate of transition from the i-th to 
the j-th odd states. Equations (§) are not independent, so they have to be supplied 
with the normalization condition P e + J2j Pj — 1- Current in the steady state equals 
to the electron flow through, say, the left barrier, and for positive V it is given by 



• r , -,)/' • ^E l ?V , r (9) 



Transition from the j-th odd state into the even state occurs when \xl > U + Ej. 
The amplitude of this transition is calculated by first order perturbation theory in the 
tunneling Hamiltonian (Q). Fermi's golden rule yields 

r«l e = 9L ^e(» L -e 3 -u), (io) 

where gi is the dimensionless conductance of the left tunnel barrier per one spin, Uj = 
(l+£j/ej)/2 is the coherence factor, 6{x) is the unit step function, and pLj = Q\ipj(r L )\ 2 , 
where Q is the volume of the grain and V'i( r i) is the value of j-th single-particle wave 
function at the left point contact r^. Energies e, are given by Eq. (|7|) and U = U + is 
defined in Eq. (|j). Similarly, the rate of transition from even state to i-th odd state, 
by tunneling of an electron from the dot to the right lead, is given by 

^ = 0*^0(^ + 6,-/*), (11) 
znti 
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where gn is the dimensionless conductance of the right tunnel barrier, Vi = (1 — £j/ej)/2, 
and pRi = Q\ipi(r^)\ 2 , where tr is the position of the right point contact. 

A change in the occupation configuration of the odd states occurs via inelastic 
cotunneling 13 . This mechanism is a virtual process in which an electron tunnels into 
j'-th available level and another electron tunnels out from the z-th level. Calculating 
this rate by second order perturbation theory in the tunneling Hamiltonian, one obtains 

r = 9L9Rd 2 u)vlp Li p m {eV - gj + Sj) 

8w 3 h(U + E j -p L )(U + E l -p R ) 1 ' 

for eV > Sj — Ei, pl < U + £j, Pr < U + £j, and zero otherwise, diverges in the 

limits pl —>■ U + Sj and p,R —> U + E{. It signals that a real transition takes over the 
virtual one. The region of applicability of Eq. fll2|) is, therefore, U + Sj — pl > 7 and 
U + Si — Pr > 7 where 7 ~ gd/4n is the width of a single-particle level in the dot 
due to the coupling to the leads, g = gi + 9r- However, the interval of biases where 
Eq. ( |T2D is not valid is narrow, and to the leading approximation in our results will 
be independent of this broadening. 

Let us now apply Eqs. (§) and (|9]) to describe the appearance of the low-voltage 
substructure of the first peak. We will consider two situations: (i) small voltage such 
that only one subresonance can emerge on the shoulder of the leading one, and (ii) 
large voltage, d <C eV < A, where the substructure of the main resonance consists of a 
large number of subresonances. 

In the first case, the chemical potentials of the left and right leads are such that 
transport through the grain involves only two levels: Eq and E\ < 6q corresponding 
to the ground and the first excited states of the odd grain. We solve Eqs. (Bp for 
probabilities Pq,Pi and P e using Eqs. (P HT^) . There are two distinct regimes of the 
source-drain voltage: (1) pi < U+eo where transport is dominated by cotunneling, and 
(2) Pl > U + Eo where state "0" is available for resonant tunneling. The substructure of 
the first resonance in the differential conductance appears in the first regime. Below we 
show that as pl passes through U+Ei, see Fig. 7, there is a discontinuity in the current- 
voltage curve. In the first regime, the total current to the leading approximation in 
9 Li gR is a sum of two contributions, I ~ I eq + I ne . The first, I eq = er ^o+2er ^i, is the 
equilibrium current coming from cotunneling. The second contribution is associated 
with the nonequilibrium population of state "1" and is given by 

C ri^o-rw+2ri^i-2r ^o ,, r <f/_|_ £l 
I ne = 2eT ^x\ r« 



2 v 1 + ?S!j Pl>U + £i 

Assuming that the voltage drop eV = Pl — Pr is larger than the energy difference 
d = Eq — Ei, the jump in the nonequilibrium current is: 

5I ne = Cl e^£ ( 1 - 4?) eV~2d (13) 
8n 2 hd \ eV) V ; 

where c\ = ^u^v\plqPri is a constant of order unity. This jump in the nonequilibrium 
current leads to the peak in the differential conductance spectra. Formula flT3"| ) has 
simple interpretation. Up to numerical prefactors it is a product of two factors: first 
is the probability of finding the grain with an unpaired electron in state "1". It is 
proportional to gn(d/d)(l — d/eV), and increases with the voltage V and as d = Eq — 
Ei — > 0. The second factor is associated with the rate in which the state "1" is filled 
with an electron, egLd/h. 
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Figure 7. Inelastic cotunneling process can drive an "odd" superconducting grain out of its ground 
state. In the ground state, the single-particle level indicated by U + £o is occupied by one electron. 
Excited states are those in which the unpaired electron is shifted to other single-particle levels. In a 
noncquilibrium steady state, low single-particle levels become available for resonant tunneling, 
leading to a subresonances structure of the differential conductance shown in Fig. 6. State j shown 
to be filled with two electrons should be understood as a coherent superposition of double occupied 
and empty states with weights i/| and respectively. 



The magnitude of the jump flT3| ) should be compared to the jump in the current 
as hl increases above U + €q, and real transition via the even state become allowed. To 
the leading order in gi and gn, the current in this regime is 

1 = ° 2e T? Ta — V' > U + e , 14 

where C2 is a constant of order unity having structure similar to C\. Comparing the 
current jump, 5I ne , with that associated with the resonant tunneling, 51, we find 

Sine 9L + ^9R t/ T 

01 OTl z 

Thus nonequilibrium population of the excited level of the odd-grain leads to the ap- 
pearance of a subresonance at small V, however, its height is much smaller than that 
of the main resonance. 

We turn now to the second regime of the parameters, d eV < A, in which many 
levels contribute to the transport. Again, we focus our attention on the cotunneling 
regime, /i£ < U + e . We show that the characteristic amplitude of the subresonances 
in this regime may become comparable to the amplitude of the main peak. 

To the leading order in gL,gn, and d/A, the steady state solution of the rate 
equations at [Ml — U + e\ + is Pq ^ 1, while for the other probabilities we have 

n ~ r (o) ' ^~ A r U) +A r U) n ' 

1 e— >o J- o—*e J- o— »e 



The characteristic number of states contributing to the current @ is large as a/ AeV jd 
so that mesoscopic fluctuations of the tunneling rates and of the inter-level spacings 
may be neglected. Additional large factor, a/ AeV / d, comes from the summation over 
the levels in Eq. flT5|), and we find 

j _ e 2 g L g R VA ( d < eV < A 

2n 2 h eo-Ej 7 { li L = U + Ej + ' 
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Once again, the current jumps each time hl passes through U + ej. This jump for large 
j (but still such that U + Sj — Ht <C eV) scales as 1/j 3 , and the ratio of the jump at 
j — 1 to the jump at the resonance level flT3| ) is given by 



51, 



51 



ne 



(g L + 4g R ) eVA 
8vr 2 dd 



d < eV < A. 



Noticing that d = Eq — S\ ~ d 2 /2A, we see that the first subresonance becomes compa- 
rable in height to the main one at voltages as small as eV ~ 4ir 2 d 3 /A 2 (gL + Ag^). 

We conclude by comparing the above results with the experimental data of Ref . [6] . 
There A « 5d, the conductances in the normal state are gn ~ lOgi ~ 1/8, and the 
leads are also superconducting. The singularity in the density of states of the leads 
imply that the effective conductance is increased by factor of 2 — 3. Neglecting inelastic 
processes and the Josephson coupling, these parameters imply that when eV ~ 2d the 
ratio of the subresonances amplitude to that of the first resonance is of order one, while 
at eV ~ 2c? « d/5 it is of order of 1%. It implies that first subresonance peak associated 
with tunneling into state "1" cannot be resolved, both because its amplitude and its 
distance from the main peak are too small. However next subresonance appear already 
at distance of order d from the main resonance, and for V > 2d have an amplitude 
comparable with the main resonance. 

4 Relaxation processes 

Central to our analysis is the assumption that the steady-state occupation config- 
urations of the electrons in the dot are far from equilibrium. This condition holds when 
the rate 1 / Ti n of inelastic relaxation processes is smaller than the tunneling rate of an 
electron into and out of the dot, l/r tun . In the opposite limit, l/r in > l/r tun , the system 
relaxes to equilibrium between tunneling events, and the electrons effectively occupy 
only one configuration. In this case one expects each resonance cluster to collapse to 
a single peak. This behavior is illustrated by the dashed line in Fig. 5 where a large 
inelastic relaxation rate \jr in = 5/r tun was included in the detailed-balance equations. 

The data shown in Figs. 2 and 3 indicate that the metal particle in the experi- 
mental system is indeed in a strongly nonequilibrium state. It is useful, however, to 
consider the various relaxation processes in our system in order to delimit the expected 
nonequilibrium regime. Relaxation of excited Hartree-Fock states may occur due to: 
(1) electron-electron interaction in the dot beyond Hartree-Fock; (2) electron-phonon 
interaction; (3) Auger processes in which an electron in the dot relaxes while another 
one in the lead is excited; (4) relaxation of an electron in the dot as another electron 
tunnels out to the lead; (5) thermalization with the leads via tunneling. The last two 
processes are small corrections since they clearly happen on time scales larger than the 
tunneling time. 

In Ref. [16] it was shown that excited many-body states of closed systems with 
energy e smaller than (g/ log g) l ^ 2 d are merely slightly perturbed Hartree-Fock states. 
In other words, the overlap between the true many-body state and the corresponding 
Hartree-Fock approximation is very close to unity. This justifies the use of our model for 
the low energy resonances since g ~ 5 therefore the energy interval 0<e< (g/ log gY^ 2 d 
contains at least the first few excited states. At high source-drain voltage, however, 
when the dot is excited to energy g l l 2 d < e < gd, tunneling takes place into quasi- 
particle states of width e 2 /(g 2 d) 17 . This width is larger than the typical separation 
between nearby resonances but smaller than d. Therefore, electron-electron scattering 
will obliterate the fine structure of resonances for high energy excitations of the dot. 
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Consider now the electron-phonon interaction. The temperature, 30 mK, is much 
smaller than the mean level spacing, therefore, the probability of phonon absorption 
is negligible, and only emission may take place. The sound velocity in aluminum is 
v s = 6420 m/sec, therefore the wavelength of a phonon associated with relaxation of 
energy uo ~ d = 1 meV is approximately 50 A, the same as the system size. In this 
regime, we estimate the phonon emission rate to be 

1 /2 \ 2 u 3 Td 

T e -ph V3 / 2ph 4 V 5 s ' 

where ep is the Fermi energy (11.7 eV in Al), and p is the ion mass density (2.7 g/cm 3 
in Al). This rate is that of a clean metal but reduced by a factor of rd/h where r 
is the elastic mean free time 18 . In ballistic systems, r is the traversal time across the 
system of an electron at the Fermi level. Assuming ballistic motion this factor is of 
order 10~ 3 . The resulting relaxation rate for u = d is therefore of order l/r e _ p /j as 10 8 
sec -1 which is similar to the tunneling rate l/r tun 6 • 10 s sec^ 1 (corresponding to 
a current of 10~ 10 A through the particle). Thus, by increasing the resistance of the 
tunnel junctions one should be able to cross over to the near-equilibrium regime shown 
by the dashed line in Fig. 5. 

Relaxation due to Auger process is estimated to be negligible. Two factors reduce 
this rate considerably: (1) it is exponentially small in w/x where w is the width of 
the tunnel junction and x is the screening length; (2) interaction between electrons on 
both sides of the tunnel junction can take place only within a very limited volume. 

5 Conclusions 

In this review it was shown that the low-voltage tunneling-resonance spectra of a 
ultrasmall metallic grains, normal as well as superconducting, reflect nonequilibrium 
electron configurations. These configurations are reached by resonant tunneling as 
well as inelastic cotunneling. The first tunneling resonance develops a substructure on 
energy scales of order of the single-particle mean level spacing, d, while high resonances 
split due to electron-electron interactions and appear in clusters of width d/g. The 
latter phenomenon is a result of electron-electron interaction beyond the orthodox 
model 1 . Relaxation due to electron-phonon interaction, which becomes important for 
high resistance tunnel barriers, will collapse the clusters. This effect can be used to 
probe the electron-phonon relaxation rate in nanometer size metal particles. 
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Appendix A 

The purpose of this appendix is to calculate the second moment of off-diagonal 
matrix elements of the interaction potential U(r — r'), and show that UijM are small as 
1/g, where g is the dimensionless conductance. The subject was discussed in several 
papers, 19 ' 20, 21 and is presented here for completeness. 
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U m = — / rfr^(r)^*(r)^(r)^(r), (16) 



When calculating off-diagonal matrix elements of the interaction potential, it 
is important take into account screening effects. The relevant two-particle interac- 
tion potential is not the bare one, t/(k), but rather the statically screened potential: 
U s (k) = f/(k)/[l + 2i/£7(k] where U(\t) is the Fourier transform of the bare two-particle 
interaction U(r — r'), and v is the density of states per unit volume. The contribution 
to the off-diagonal matrix elements comes only from spatial fluctuations in the electron 
density (non-zero modes) for which screening is established at very short time, of order 
of the time it takes for a plasmon to propagate through the system. The latter is much 
shorter than the relaxation time, t c , of fluctuations in the electron density. Thus for 
large v the screened interaction potential, U s (r — r'), is close to a 5-function. 

Consider, therefore, the integral 

1 

where no two indices are the same. Clearly on average (Uijki) = since wave func- 
tions associated with different eigenvalues are independent and (ip) = 0. To estimate 
the magnitude of the off-diagonal matrix elements we calculate the second moment 
(\Uijki\ 2 )- The square of the matrix element, \Uijki\ 2 , contains four pairs of wave func- 
tions in the form ip*{r)ipi{r'). Since the correlation between wave functions and eigenen- 
ergies are only to order 1/g, one can approximate these pairs as 

V*(r)Vi(r') « ^~ [G(r, r'; E, - ir,) - G(r, r'; E, + irf)} , (17) 

where d is the single-particle mean level spacing, rj is a positive energy which will be 
taken to zero at the end of the calculation, and G(r, r'; E) is the single-particle Green 
function at energy E. Two basic correlators emerge when calculating the ensemble or 
the energy average of \Uijki\ 2 - These correlators, known in disordered diagrammatic 
nomenclature as the diffuson and Cooperon, are 

Il^r, r') = (G(r, r';E + LU + ir])G(r\ r; E - ir])), 

and (G(r, r'; E + ir])G(r', r; E + u — ir])). For systems with time reversal symmetry, 
considered here, these correlators are the same. In the semiclassical limit, 



II w (r,r , ) = 27ri/X; 



where the sum is over all classical relaxation modes, i.e. diffusion modes in the case of 
disordered grains and Perron- Frobenius modes in chaotic systems^ 7 M are the corre- 
sponding eigenvalues, and x M (r) [xAt( r 0] i s the projection of the Perron- Frobenius left 
[right] eigenfunctions on the real coordinate space at fixed energy E. 

With the help of fllED, (|T7|) and (|IBD, and assuming all energy differences (such as 
Ei — Ej) to be much smaller than h'ji, one obtains 

/,/\-' 

\Uijki\ ) 

where g= hji/d is the dimensionless conductance of the system, and d is a constant of 
order unity given by 




7m/7i Ivhi 



3 For simplicity we consider here chaotic systems in the form of billiards, namely the Hamiltonian 
contains only a kinetic part, and chaotic dynamics is due to the irregular boundary. 
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Here we assumed for simplicity that 71 is real. In case it contains also an imaginary 
part, the same formula applies with the substitution 71 — > Re{7i}. Notice that there 
is no zero mode contribution to (\Uijki\ 2 ), since only density fluctuations associated 
with non zero- modes can induce scattering and contribute to Uijki- Mathematically 
this results from the fact that eigenfunctions, Xo( r ) an d Xo( r ) ; associated with the zero 
mode are real, and since 11^ is always calculated at a finite energy deference, u>, taking 
its real part excludes the zero-mode contribution. 

The rest of this appendix is a semiclassical derivation of formula (0). We begin 
by writing Green's function in the semiclassical approximation as a sum of two terms 22 : 



G(r, r'; E ± i V ) ~ G (r, r'; E ± i V ) + y ^ £ A^^^'.^l 

Here Go is the Weyl contribution associated with "zero length" trajectories. This 
term is important only at distances |r — r'| of order of the particle wavelength and 
therefore can be neglected. The second term is a sum over all classical trajectories 
from r' to r, in which / is the number of degrees of freedom, S rr > t i(E) is the action, 
and A rr / t i is the corresponding probability amplitude. It is convenient to introduce a 
local coordinate system in which r is the time along the trajectory, and r_i_ are the 
coordinates perpendicular to the trajectory. In these coordinates 22 



where f and r' denote the velocity of the particle at the final and initial points, and 
p' ± is the conjugate momenta to r'±. 

Expanding S rr ' t i(E + uj) to first order in lu, and using the diagonal approximation 
for the product of the two Green functions, one obtains 

iWr, O = ^ E lArr'A'e^ = ^ / dtP(t)e^\ (19) 
where 7] = dS rr ^i{E)/dE is the time associated with the Z-th trajectory, and 



p<«> = E>t- l (fe) ( '<«-i». 



(20) 



Next we show that P{t) is the projection of the classical propagator in phase space 
onto configuration space at fixed energy E, namely 

P{t) = J dp' Jdp5[E-H(r,p)\ (r,p|e- |r',p'), (21) 

where H(r, p) is the classical Hamiltonian of the system and e~ ct is the evolution 
(Perron-Frobenius) operator for time t. In the coordinate system introduced above, 
the Hamiltonian function H is the conjugate momentum to the time coordinate r 
along the trajectory, therefore 

f d M f d J-f 1 

Pit) = j — d P± J —dpfJ(E - H)8(H - H' t )8(r - r^5(r ± - r\ t )5(p ± - p' ±t ) 
-^^(t-Ti)j v dp' x S(v ± -v' ±t ), 
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where subscript t denotes the value of the corresponding coordinate after time t starting 
from the phase space point (r', p'). Since the energy of the particle is fixed, the integral 
reduces to a discrete sum over trajectories from r' to r. The contribution to each 
trajectory, I, comes from an infinitesimally small region of the coordinate p' ± denoted 
by Ti. Straightforward integration yields the result (^Hp. 

Starting now form ([H]) and using the spectral decomposition of the Perron- Frobenius 
operator 

(r,p|e- £ y,p'> = f5(pf -p'0£e-^(r,n)^(r',n'), 

where n denote the direction of the momentum, while (p^ and ip^ are the left and right 
eigenfunctions, we get 

P(t) = ^£e->* f ^^(r,n)^(r',n') = h f vJ2^%(*)x»(r')- 

Here Qf — J dn is the solid angle of a sphere in / dimensions, and v — J dp5(E — H)/hf 
defines the density of states per unit volume. Substituting this expression in (|T9|) and 
performing the time integration yields the required result (fjjj). 



Appendix B 

In this appendix we derive formula (|7|) for the cost of introducing an additional 
electron into an odd state: Ej = E 2m — -Esm-i (ignoring the charging energy). Denote 
by Fn = E^- — [inN the free energy of the system where En is the ground state energy 
with N electrons, and /xjv is the chemical potential. For the superconducting state with 
unpaired electron in level j one has 



= £(& - e*) + A 2 /A + 0, e k = ^ + A 2 , (22) 

where A is the pairing coupling constant. In the intermediate state where another 
electron tunneled into the j-th state and pairs with the already existing one: F 2m = 
J2k(£k ~ e fc) + A 2 /A. The parameters in these two equations, A and \i (which are 
functions of N), are determined by the relations: 

dF N BF N 



dA N d/j, 



N 



Notice that the single-particle energies are measured with respect to the chemical 
potential, thus for the differentiation with respect to /i^ it it is convenient to introduce 



£fc = £fc - Vn, and efc = V (Cfc - Vn) 2 + Afy, where £ fc are independent of fi N . 
The above derivatives for the even case, N = 2m, give 

E / ~ 1 =T» 2m = ^fl- . _ ~ ik ~^ m 1, (23) 



k 2J^ k -/2 2m ) 2 + A 2 2m X k V J(Zk-fi2m) 2 + & 



2m, 



while for the odd case, N = 2m — 1, 

E 1 , .... _ _ . 

(24) 



1 n o I i P2m-1 
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We expand e» 



E% m — E. 



2m- 1 



F2m{^2m\ A 



2m J 



^2m-l(^2m-l) ^2m-l) + 2m// 2m — 



(2m — l)/x 2 m-i t° linear order in the differences A 2m — A 2m _i, and // 2m — /i 2m _i: 

£j ~ F 2m (fi2m] A 2m ) — F 2m -i(fi2m] A 2m ) + /i 2m 



/4m- 1 V(£. 



A*2r, 



) 2 + A 2 



2m ' 



(25) 



We are left, now, with the problem of finding the dependence of /4m-i on the unpaired 
electron state j. For this purpose we choose // 2m as the reference point and calculate 
the difference /4m- i 



2m -2 




//2m- From Eqs. (P3|) and (gj) we have 

4 



..0) 



£fc — //2m 

(6c - At2m) 2 + A 
/ 



ij — //2r, 



2 

2m / 



A 2 

^2m, 



— (/4m- 1 — A t 2m) X! 



-1 



) 2 + A 2 



+ 



2m 



((a.-/^2m) 2 + A 2 m ) 



\ 3/2 



(The correction to A vanishes in the limit where the single-particle mean level spacing 
d is much smaller than the superconducting gap, (i< A). Since the first sum in the 
second line equals 2m we have 



(£j ~ l^2r, 



) 2 + A 2 



+ E 



(/4m-l ~ / i 2m)A 



2m 



2 rn 



A 



+ (/4m- 1 A*2ri 



^-/^2m) 2 + A 2 m j 

d£ A 2 
T(^2 + A 2)3/2' 



\ 3/2 



Thus 



A*2r, 



(?) ^ I a 

fhm-i + 2 1 A 



and substituting this result in fl25|) we obtain (|7]). 
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